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Abstract. We investigate the dependence of the order parameter profile, local and 
total susceptibilities on both the temperature and external magnetic field within the 
mean-filed Ginzburg-Landau Ising type model. We study the case of a film geometry 
when the boundaries of the film exhibit strong adsorption to one of the phases 
(components) of the system. We do that using general scaling arguments and deriving 
exact analytical results for the corresponding scaling functions of these quantities. 
In addition, we examine their behavior in the capillary condensation regime. Based 
on the derived exact analytical expressions we obtained an unexpected result - the 
existence of a region in the phase transitions line where the system jumps below its 
bulk critical temperature from a less dense gas to a more dense gas before switching on 
continuously into the usual jump from gas to liquid state in the middle of the system. It 
is also demonstrated that on the capillary condensation line one of the coexisting local 
susceptibility profiles is with one maximum, whereas the other one is with two local 
maxima centered, approximately, around the two gas-liquid interfaces in the system. 

Keywords: Exact results, Phase transitions and critical phenomena. Classical phase 
transitions (Theory), Finite-size scaling. Phase diagrams (Theory) 


PACS numbers: 64.60.-i, 64.60.De, 64.60.F-, 75.40.Cx 



Exact results for the temperature-field behavior of a mean-field model 

Contents 


2 


1 Introduction 3 

2 The Ginzburg-Landau mean-field model and basic expressions defining 

the finite-size quantities of the system 5 

3 Analytical results for the scaling behavior of the order parameter 

profiles 7 

3.1 Analytical representation of the order parameter profiles in the case of 

zero field. 7 

3.2 Analytical results for the order parameter profiles in the case of nonzero 

field . 9 

4 Analytical results for the scaling behavior of the susceptibility 13 

4.1 Exact results for the local susceptibility profiles in the case of nonzero field 13 

4.2 Exact results for the local susceptibility profiles in the case of zero field . 15 

4.3 Exact results for the total susceptibility in the case of nonzero field ... 16 

4.4 Exact results for the total susceptibility in the case of zero field. 18 

5 Discussion and concluding remarks 18 

Appendix 19 


References 


21 






Exact results for the temperature-field behavior of a mean-field model 

1 . Introduction 


3 


The understanding of the phase behavior of fluids confined in narrow regions, including 
the fluid mediated interactions, is of crucial importance in the physics of fluids in porous 
media, for colloidal physics, for many applications and modern technologies such as 
lubrication, adhesion or friction, in micro and nano-fluidics as well as for the proper 
interpretation of surface force experiments 

In the current article the order parameter prohle and its response functions to 
an externally applied ordering held h, as well as the free energy will be investigated 
as fnnctions of both the temperature T and h for the three-dimensional continnnm 
mean-held Ising model with a him geometry cx)^ x L. We will consider the full (T, h) 
plane for the case when the bounding surfaces of the system strongly prefer the ordered 
phase of the system which can be a simple hnid, magnetic system close to their 
respective critical point, or a binary liquid mixture close to its demixing point. This is a 
standard model within which one studies phenomena like critical adsorption [1, 12-24], 
wetting or drying [21, 22, 25-28], surface phenomena [29, 30], capillary condensation 
[1, 16, 17, 19, 26, 31, 32], localization-delocalization phase transition [33-35], hnite-size 
behavior of thin hlms [2, 16, 33, 35-42], the thermodynamic Casimir ehect [19, 43-45], 
etc. One normally derives the resnlts for h = 0 analytically [43, 45, 46] while the h- 
dependence is studied numerically either at the bulk critical point of the system T = Tc, 
or along some specihc isotherms - see, e.g., [19, 24, 34, 44, 46-48]. 

In the current stndy we will present analytical results for the h-dependence of the 
model and will provide results for the full (T, h) dependences of the order parameter, 
local and total susceptibilities in the (T, h) plane. In what follows we will mainly use 
the magnetic terminology but when considering capillary condensation we will also use 
that one of a simple fluid system in order better to reflect the physics of the obtained 
exact mathematical results. 

In order to be more specihc let us remind some facts and dehnitions pertinent to 
the above mentioned problems. 

If a huid, or a magnetic system possesses a bounding surface its phase behavior 
as a function of its temperatnre T, excess chemical potential A/i oc h and the material 
characteristics of the surface are essentially enriched - near the snrface one can have, 
e.g., phenomena of wetting or drying [27]. In the vicinity of the bnlk critical temperature 
Tc of the bulk system, one observes a diversity of surface phase transitions [29, 30] of 
different kind in which the surface orders before, together, or after ordering in the bulk 
of the system, which is known as normal (or extraordinary), surface-bulk and ordinary 
snrface phase transitions. For a simple hnid or for binary liquid mixtures the wall 
generically prefers one of the huid phases or one of the components. In the vicinity of 
the bulk critical point the last leads to the phenomenon of critical adsorption [1, 13-24]. 
Obviously, the surface breaks the spatial symmetry of the bulk system. The penetration 
depth of the ehects due to the existence of a bounding surface in the body of the system is 
set by the correlation length ^ of the order parameter; ^ becomes large, and theoretically 
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diverges, in the vicinity of the bulk critical point (Tc, h = 0): f(T —)■ Tfi, h = 0) ~ 

T = {T — Tc)/Tc, and f{T = h —)■ 0) ~ fo^h\h/{kBTc)\~‘'^^, where u and A are the 
usual critical exponents and and fo,h are the corresponding nonuniversal amplitudes 
of the correlation length along the r and h axes. When at least one of the spacial 
extensions of the system is hnite, as in the cx)^ x L him geometry we consider in the 
current article, one terms the corresponding system a hnite system. If in such a system ^ 
becomes comparable to L, the thermodynamic functions describing its behavior depend 
on the ratio L/^ and take a scaling form given by the hnite-size scaling theory [2, 33, 40- 
42]. One observes, inter alia, shift of the critical point of the system [16, 35, 37-39] from 
Tc to Tc^l- When in the hnite system there is a phase transition of its own Tc^l is a 
true critical point. Below Tc, if the conhning walls of the him are of the same material, 
the capillary condensation occurs [1, 16, 17, 19, 31, 32] where, e.g., the liquid vapor 
coexistence line shifts away from the bulk coexistence into the one-phase regime. This 
coexistence line of hrst order transitions ends at a point Tc^Li which is normally identihed 
with the capillary condensation point Tcap that, on its turn, is considered to correspond 
to the highest temperature at which the entire capillary hlls with liquid. In the current 
article we will demonstrate that T^ap = Tc^l is not always true. We will show that, 
within the studied model, Tcap differs from Tc^l on a scale determined by . 

As already stated above, near a conhning wall the symmetry between the two 
phases of a simple huid or between the two components of the binary liquid mixture is 
violated in that one of these phases or components is preferred by the boundary. Thus, 
the order parameter prohle, the density or the composition, becomes a function of the 
perpendicular coordinate x. This can be modeled by considering local surface helds hi 
and /i 2 acting solely on the surfaces of the system. When the system undergoes a phase 
transition in its bulk in the presence of such surface ordering helds one speaks about the 
” normal transition [49]. It has been shown that it is equivalent, as far as the leading 
critical behavior is concerned, to the ” extraordinary' transition [30, 49] which is achieved 
by enhancing the surface couplings stronger than the bulk couplings. In the remainder 
of this article we will use the surface held picture. It has been demonstrated [30] that 
when hih 2 Q for the leading critical behavior of the system is sufficient to investigate 
the limits hi, /i 2 —)■ ±oo. Obviously, there are two principal sub-cases hi = ^2 —)■ -|-oo, 
and hi = —h 2 —)■ +C )0 corresponding to hih 2 > 0 and hih 2 < 0. One usually refers to 
the former case as the (-I-, -|-) boundary conditions and to the latter case as the (-I-, —) 
boundary conditions. In the current article we will be only dealing with the behavior of 
the system under (-f, -|-) boundary conditions. For such a system the hnite-size scaling 
theory [2, 40, 50] predicts: 

• For the magnetization (order parameter) prohle 

m{z,T,h,L) = -- ahT~^l''Xra{zlL,Xt,Xh) (1.1) 


where 


xt = atrL^!'', 


Xh = ahhL^^’'. 


( 1 . 2 ) 
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• For the local (layer) susceptibility profile 

ksTxi{z,T,h,L) = ■^m{z,T,h,L) ~ {z/L,xt,Xh) (1.3) 

with X^{z/L,Xr,Xh) = ^Xm{z/L,Xr,Xh)] 

• For the total susceptibility x(T, h, L) = L~^ Xi{z, T, h, L)dz one has 

ksT x{T,h,L) zz alL^/^X {xt,xn). (1.4) 

In Eqs. (1.1) - (1.4), fi and 7 are the critical exponents for the order parameter and the 
susceptibility (compressibility), the quantities at and an are nonuniversal metric factors 
that can be hxed, for a given system, by taking them to be, e.g., at = 1 / and 

0/1 = 1/ Since the Ising system with a film geometry (xfi x L possesses a critical 

point Tc^l of its own with coordinates {x[^\x^^^) the scaling functions Xm,X^ and X 
will exhibit singularities near this point. For example 

X{xt,x^f;^) ~ Xc,t (xt - ^ , Xt^ xl''\ (1.5) 

where the subscript in 72 reminds that 72 is the critical exponent of the two-dimensional 
inhnite system that is to be distinguished from the corresponding exponent 7 for the 
three dimensional bulk system. 

In the current article we will derive new exact analytical results for the scaling 
functions Xm,X^ and X for the Ginzburg-Landau Ising type mean-held model. Let us 
recall that in the mean-held approximation /3 = z/ = l/2,A = 3/2 and 7 = 72 = 1. For 
the version of the model considered here = 1 and fo,h = 1/^ [19, 44, 51]. 

We begin our study by presenting a short dehnition of the model that will help us 
to introduce some of the notations used further in the article. 

2. The Ginzburg-Landau mean-field model and basic expressions defining 
the finite-size quantities of the system 

Let us consider an Ising type critical system in a parallel plate geometry described by 
the standard 0^ Ginzburg-Landau Hamiltonian 

X[(j);T,h,L]= [ C{(l),(fi)dz, (2.6) 

Jo 

where 

C = (fi) = ^0'^ + - h(j). (2.7) 

Here: L is the him thickness, is the order parameter at the perpendicular position 
{0 < z < L), T = (T — Tfij/Tc is the bare reduced temperature with r = 0 dehning the 
bulk critical temperature, h is the external ordering held, g is the bare coupling constant 
and the primes indicate diherentiation with respect to the variable ; 2 . Normally, one 
adds to the right-hand side of Eq. (2.6) a surface-type term with parameters used to 
impose the boundary conditions on the system. We will do that by simply requiring the 
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behavior of the order parameter at the boundaries of the system at z = 0 and z = L to 
be of a given prescribed type. 

The extrema of the functional iF are determined by the solutions of the 
corresponding Euler-Lagrange equation 


ddC_dC_ 
dz d(j)' d(j) 

which, on account of Eq. (2.7), reads 


( 2 . 8 ) 


0 " — 0 T -\- g (ffi + h = 0. 


(2.9) 


Multiplying the above equation by 0' and integrating over one obtains the following 
Erst integral of the system 

^ 0 '^ - ir 0 ^ - ^^ 0 ^ + h 0 = c, ( 2 . 10 ) 

where c is the constant of integration. 

In the present article we choose the so-called (-I-, -1-) bonndary conditions: 
lim 0 (; 2 ) = hm 0 ( 2 ;) 0 ^^ = -|-cx). Due to the symmetry, under snch boundary 

conditions one shall have that 0'(L/2) = 0. It is easy to determine the general behavior 
of (j){z) near the boundaries for any fixed finite values of r and h. Since (j){z) —)■ cx) near 
the boundaries then, say, for the left bonndary, Eq. (2.10) becomes 





Solving this equation leads to 


0(2:) ~ 


1 

- Zo 


( 2 . 11 ) 


( 2 . 12 ) 


where zo is the position of the bonndary (in the case considered zq = 0). Note that 
this leading behavior of the order parameter prohle near the boundary does not depend 
neither on r, nor on h. From Eq. (2.12) it follows that the integral in Eq. (2.6) 
diverges. Thus, when a comparison of the valnes of iF for different states of the system 
is needed either some cnt-off of the system near the bonndaries is necessary, or some 
special procednre shall be devised. 

From Eq. (2.9), using the definition for the local layered susceptibility 


Xi{z\t, h) 


(90(2:|r, h) 

mi 


one obtains that xi fulfills the relation 


( 2 . 13 ) 


-X"+ a + 3g'p)xi = 1. 


( 2 . 14 ) 


On general grounds, as well as from Eq. (2.12), for a system with (-I-,-I-) boundary 
conditions one immediately obtains xi{.z —)■ zfir^h) = 0. Combining Eqs. (2.12) and 
(2.14) leads to Xi{.z\t, h) ^ {z — Zofi/i when — )■ zq. 
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3. Analytical results for the scaling behavior of the order parameter profiles 


In terms of the scaling variables 

C = z/L, xt = TL^Ifi Xh = ^ghL^lfi (3.15) 

(j){z) = ^J^L~f^/''Xm{C,\xt,Xh), (3.16) 

with fi = V = 1/2 and A = 3/2, Eq. (2.9) for the order parameter prohle takes the form 
AKC) = A„(C) [x, + 2At(C)] - y. (3.17) 

The solutions of Eq. (3.17) determine the extrema of the energy functional 

T= [' f{XUO,X^^{0)dC, (3.18) 

Jo 

where 

/(A„(C), a;(c)) = K,(C)]" + XiK) + x,XHC) - xMC) ( 3 . 19 ) 


is the energy density. Hereafter, the primes indicate differentiation with respect to the 
variable ( which, as follows from Eq. (3.15), varies in the closed interval [0,1]. According 
to Eq. (2.10) the hrst integral of Eq. (3.17) reads 

lAKOl" = P[XA, P\Xn,] = X‘JC) + x,X^(() - x.XUO + <r.(3.20) 
where e denotes the respective constant of integration. 


3.1. Analytical representation of the order parameter profiles in the case of zero field 

When h = 0 the magnetization prohle is known exactly [43] in terms of two mutually 
related implicit equations: 

a) when Xt = tLA > — 

cln|2A-(l;)C;*;] 


.^(Cki.O) - 

where > 0 is to be determined from 
Xt = [2K{k)]\2l^ - 1); 

b) when xt < — 

* snl2A'(i)C;fc]’ 

where > 0 is to be determined from 
Xt = -[2KCk)]\P + l). 


(3.21) 


(3.22) 


(3.23) 


(3.24) 
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Here K{k) is the complete elliptic integral of the first kind, dn(C;/c) and sn{(]k) are 
the Jacobian delta amplitnde and the sine amplitnde fnnctions, respectively. The bnlk 
critical point T = Tc corresponds to k'^ = 1/2. Note, however, that within the mean- 
field theory the magnitnde of the variable 0 is not nniversal, in that it is mnltiplied 
by the nonnniversal factor \/2/g - see Eq. (3.16). Finally, we stress that the choice 
of two parameterizations (see Eqs. (3.22) and (3.24)) of the scaling fnnctions in Eqs. 
(3.21) and (3.23), is just for convenience; it allows one to avoid using imaginary values 
of k and k. Indeed, one can transfer any of the set of equations into the other one. 
For example, defining k as k = ik/k' where k''^ = 1 — k^, and taking into account the 
following properties of the elliptic functions [52, 53] K{k) = k' K{k) and 


dn(M; ik) 


y/l + k‘^ 


(3.25) 


sn{u]ik sn{uy/l+~k^-,k/y/l+~k^)' 

one can easily check that the pair of equations (3.21), (3.22) is equivalent to the pair of 
equations (3.23), (3.24). 

In order to utilize the symmetry of the problem it is helpful to move the coordinate 
frame so that the origin of the ( axis is at the midpoint of the him [46]. Taking into 
account that [52] 

dn[u-\-K{k); k] k' 


STi[u -\-K{k); k] cn(-u; fc) ’ 
from Eq. (3.21) one obtains [46] 


and cn{iu;k') = l/cn{u;k) (3.26) 


X^{C\xt) = X^o cn[t2K{k)C;k'] 


(3.27) 


where ( G [—1/2,1/2] and 

Xmo = 2k'K{k) = 2K{k). 


(3.28) 


Since cn(0; fc) = 1, one has Xmo = Xm{0\xt). Eq. (3.27) is the simplest representation 
of the order parameter prohle in a system with strongly adsorbing boundaries we are 
aware of. 



Figure 1. Color online. Plot of the scaling function of the order parameter profile 
XmiClxt) - see Eq. (3.27), normalized per its value in the middle of the system 
N,„(C = 0|xt =0) forr = Te. 
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The typical behavior of Xm{C\xt) is shown on hgure 1 for three different 
temperatures: at, well below, and well above the bulk critical temperature T^. 


3.2. Analytical results for the order parameter profiles in the case of nonzero field 


When h ^ 0 the study of the order parameter (magnetization) prohles was carried 
out numerically [44], Below, however, we shall give an analytical representation of the 
foregoing prohles by means of Weierstrass elliptic functions, which is quite similar to 
that presented above for h = 0 up to the lack of the parametrization through the elliptic 
modulus achieved in that case, see Eqs. (3.21) - (3.24). 

First, let us recall that we are interested in real-valued solutions Xm{C) of Eq. 
(3.17), corresponding to given values of the parameters Xt and Xh, which are smooth 
in the open interval (0,1) and satisfy the (-f, -1-) boundary conditions. Evidently, 
for each such solution the polynomial P[Xjn] should have at least one real root 
Xmo- Otherwise the considered solution would be either strictly increasing or strictly 
decreasing, contrary to the required boundary conditions, since its derivative X!^(() 
would be either strictly positive or strictly negative as implied by the particular form of 
Eq. (3.20). Consequently, the constant of integration e corresponding to such a solution 
can be cast in the form 


^ XmO (^mO ^tXmO . 


(3.29) 


Now, given a triple of values of the parameters Xt, Xh and X^o, each real-valued solution 
Xjn{() of Eq. (3.20) can be expressed, following [54, §20.22, §21.73], in the form 


(Cl^t) ^h) Xfno) 


^ 6 X ^0 {xt + 2 X^o) - 2,Xh 

12 p(C-|;^2,^3) - (a^t+ 6 X^o) 


(3.30) 


Here, 92 , 93 ) is the Weierstrass elliptic function, 92 and 93 are the invariants of the 
polynomial P[Xm], which according to [54, §20.22, §21.73] and Eq. (3.29) read 


92 ^mO ^t^mO ^h) , 


(3.31) 


93 


1 


[27a;^ -h 2x1 + 72a;Wmo {X^o + XtXmo -Xh) ■ 


(3.32) 


Note that each function of form (3.30) has the following important properties. First, 
its graph in the (C, X^) plane is symmetric with respect to the line parallel to X^ axis 
and passing through the point (1/2, 0). This is because the Weierstrass elliptic functions 
have the property ^{-^-, 92 , 93 ) = P 92 , 93 )■ Next, 


X. 


2 \ xt, Xh, X^o 


= X 


mO, 


XL 


2 \ xt, Xh, X^iO ) 0, 


(3.33) 


since the function p — 1 / 2 ; 92 , 93 ) has a second-order pole at C = 1/2, i-e., 
limp(/' — 1/2; 92, ~ + 00 . Finally, the considered function has a local minimum 

at C = 1/2 if and only if 


2X^0 {xt + 2 X^ 0 ) -Xh>0 


(3.34) 
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since Eqs. (3.17) and (3.33) imply 

^ [2X„o {xt + 2X^o) - Xh] ■ (3.35) 

Let us stress that not any function of the form (3.30) satisfies the required (+, +) 
boundary conditions. Actually, a function of form (3.30) corresponding to a given couple 
of values of the parameters Xt and Xh meets these conditions if and only if X^o is such 
that: 

(a) the denominator of the second term in the right hand side of expression (3.30) attains 
zero at C = 0 and C = 1) which in view of the relation p (—^; g 2 , gs) = p g 2 , gs), means 

ip{xt,Xh,Xmo) = 0, (3.36) 


where 


(pi^Xf, X}i^ Xffig') . 12p 



(3.37) 


fbj Xmo satisfies constraint (3.34), meaning that the function Xm {C\xt, Xh, X^o) attains 
its minimal value at C = 1/2, i.e., at the center of the system. 

Thus, in order to determine the smooth for each ( G (0,1) functions 
{C\xt,Xh,Xmo) of form (3.30) satisfying the considered (+,+) boundary conditions 
for given values of the parameters Xt and Xh one should find all the solutions X^o of 
the respective transcendental equation (3.36), which are such that constraint (3.34) 
is fulhlled. Any such function will represent an order parameter prohle satisfying the 
(+, +) boundary conditions. In the cases in which the parameters Xt and Xh are such that 
there is more than one value of the parameter X^o satisfying the above requirements, 
i.e., there is more than one order parameter prohle satisfying the (+,+) boundary 
conditions, on physical grounds we chose the one that minimizes the truncated energy 


Str{xt,Xh,Xmo) = 2 f{X^{C),X'^{C))dC, (3.38) 

were <5 is a small positive number, i.e., 5 1. We work with the truncated, instead 

with the full energy, since the integral (3.18) determining the energy S of the states of 
the system is divergent. Below we justify this procedure in a mathematically rigorous 
way. Before proceeding to that, let us introduce an approximation of the considered 
order parameter prohles near the singular point C = 0- 

Given xt, Xh and e = —X^o + XtXmo — Xh), the solution of Eq. (3.20), which 
is unique, can be approximated near the singular point C = 0 by the function 


XmiO 


1 _ ^ . Xh.2 , 73;^ - 36£ 3 _ XfeXt 4 
c 6 ^ ^ 8 ^ 360 ^ 48 ^ 


(3.39) 


Indeed, the substitution Xm = Xm in Eq. (3.20) gives 


a:„(c) 


- x^(C) - x.xf^iC) + x.XmiC) - ^ = ( + ^ - 


ext 


3x1 


C" + O(C=).(3.40) 
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Consequently, using the approximate solution Xm of Eq. (3.20), given by Eq. (3.39), 
one can approximate the energy density (3.19) near the singular point C = 0 as 

fixuci KiO) = ^ + ^- f- ^ + | + ^(0- (3.41) 

Next, using the symmetry of the energy density with respect to the point ( = If2, 
implied by the symmetry of the general solution (3.30) with respect to this point and 
Eq. (3.19), we rewrite the energy (3.18) in the form 

£{xt,Xh,Xmo) =2 [ f {Xm{C),X'^{())dC + Strict, Xh,Xmo;S). (3.42) 
Jo 

Now, let Xmoi and Xm 02 be two different values of the parameter Xmo determining two 
different states of the system for same values of the parameters Xt and Xh- Then, using 
near the singular point C = 0 fhe approximate energy density (3.41) corresponding to 
the approximate solution (3.39), we observe that 

S{Xii Xh^ XYaof) £ X}i, X^nOl) £tr{,XtT Xfi, X^q2j S') £fj-(^Xf, Xj^, X^qi^ S') (3.43) 

2 

+ g Xfi, Xm02) — ^{xt, Xh, XmOl)] S + 0(6'^) 

and, hence, the difference between the energies of the two regarded states is well dehned 
and determined by the difference between the two respective truncated energies up to 
terms of order 0(5). As it is clear from Eq. (3.43), S determines the precision with which 
we determine the energy differences between any two solutions of the oder parameter 
problem. Since the value of S is on our disposal, we can, at least in principle, determine 
these energy differences to any prescribed precision. A numerical verihcation of the 
above relation is presented in the Appendix. 

In the remainder, using the derived exact analytical expressions described above 
in this Subsection, we study the behavior of the order parameter prohles in the critical 
and in the capillary condensation regimes. It should be stressed, however, that the 
solutions Xmo of the transcendental equation (3.36) corresponding to given values of 
the parameters Xt and Xh are determined numerically and those of them that meet the 
other necessary conditions are identihed by inspection. 

Our Erst observation concerns the number of the solutions of the considered (+, +) 
boundary value problem. There are values of the parameters Xt and Xh for which we 
find only one solution of the problem, but there are also regions in the temperature-held 
plane where there exist three solutions satisfying the (-I-, -|-) boundary conditions, which 
reduce to two in certain limiting cases. These alternatives are illustrated in figures 2 
and 3. In figure 2 we show the evolution with xt (left) and Xh (right) of the value 
Xmo of the order parameter in the middle of the him for three diherent values of the 
temperature Xt and for one value of the held Xh, respectively. The order parameter 
prohles corresponding to Xh = —250, Xt = —30 and Xh = —250, Xt = —32 are depicted 
in figure 3, where the stable prohles are represented by thick curves. 

The existence of more than one order parameter prohles for one and the same 
temperature-held combination is a necessary, but not sufficient, condition for the 
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Figure 2. The dependence of the value of the order parameter in the middle of 
the film on Xt and Xh- Left: Xt = —20 (dotted), Xt = —33.8105 (thick) and Xt = —40 
(dashed); right: Xh = —250. 




Figure 3. Order parameter profiles at Xh = —250 corresponding to Xt = —30 (left) 
and Xt = —32 (right), where the stable profile is represented by the thick curve. 


occurrence of capillary condensation transition. The latter takes place when at least 
two of the observed order parameter profiles have the same energy meaning that they 
coexist. 

Based on the derived exact analytical expressions we obtain the phase diagram (see 
figure 4, left) and an unexpected result of the existence of a curve in the temperature-field 
plane (see figure 4, right) where the system jumps below its bulk critical temperature 
from a less dense gas to a more dense gas (see figure 5) before switching on continuously 
into the usual jump from gas to liquid state in the middle of the system in the capillary 
condensation regime (see figure 6). Some technical details related to the determination 
of the data presented in these plots are given in the Appendix. 

We stress here that we term a given state of the system ”gas”, or ”gas-like”, if 
Xm(l/2\xt,Xh,Xmo) < 0 and ’’liquid”, or ”liquid-like”, when Xmil/2\xt, Xh, Xmo) > 0. 
We remind that under the (-f, -b) boundary conditions studied in the current article 
one always has Xm{C\y^t,Xh,^mo) > 0 for ^ close enough to 0 or 1, i.e., one always 
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Figure 4. Phase diagram (left) showing the border line crossing which the system 
jumps from a gas to either a denser gas or liquid state in the middle of the 
system. This curve ends at the critical point of the finite system = 

(-25.6983391,-236.0350005). The figure on the right shows the pre-capillary- 
condensation curve where above and below the system jumps from a less 
dense gas to a more dense one. For T < Tcap the order parameter in the middle of the 
system jumps from a gas to a liquid state. 


observe a ” liquid-like” state near the boundaries of the system. When one lowers the 
temperature the following is happening. Above one has a single order parameter 
prohle that satishes the (-I-, -|-) boundary conditions. Near the phase transitions line 
there are already three such prohles two of which provide at a point belonging to the 
phase line and characterized with given xt and Xh the minimum of the energy of the 
system, i.e., they describe the phase coexistence between the gas and the liquid phases. 
When crossing this line the ” liquid-like” order parameter prohle changes abruptly with 
the liquid phase intruding deeper into the capillary. It turns out, however, that in 
a given temperature range the two liquid branches stemming from the two surfaces 
of the capillary do not meet in the middle, but a ”gas-lke” gap still exists there for 
temperatures T^ap < T < T^^l- We call the phase coexistence line for this special 
case pre-capillary-condensation curve. Upon further reduction of the temperature or 
increase of the magnitude of the negative external held the density of the huid in this 
gap continuously increases reaching its liquid value. Thus, if one dehnes the capillary 
condensation temperature Tcap as the highest one at which the entire capillary hlls with 
liquid one will obtain Tcap that dihers from on a scale determined by . Of 
course, this is a result that follows within the model considered and will be desirable to 
check if it is a specihc feature of the model or if it can be experimentally verihed. 

4. Analytical results for the scaling behavior of the susceptibility 

f.l. Exact results for the local susceptibility profiles in the case of nonzero field 


With respect to local susceptibility 

Xi{z\xt,Xh) = L'f/’'X^{C\xt,Xh), 


(4.44) 
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0.0 



0.2 0.4 0.6 0.8 1.0 


0.0 0.2 0.4 0.6 0.8 1.0 


Figure 5. The pair of order parameter profiles coexisting at pre-capillary-condensation 
for Xh = —236.2 and Xt = —25.795. 



0.0 0.2 0.4 0.6 0.8 1.0 



0.0 0.2 0.4 0.6 0.8 1.0 


Figure 6. The pair of order parameter profiles coexisting at capillary condensation 
for Xh = —250 and Xt = —33.8105. 


with 1/ = 1/2, 7 = 1 in the model under consideration, from Eqs. (2.13) and (4.44) one 
derives 


d 


(4.45) 


In particular, for the local susceptibility in the middle of the system one has 

X^{l/2\xt,Xh) = 2Xmo{xh), (4.46) 

where the dot indicates differentiation with respect to the variable Xh- According to Eq. 
(2.14), X-^ satishes the equation 

A'"(C) - [x, + 6A7(C)] .Aj(C) = -1. (4.47) 

On the other hand, differentiating the hrst integral (3.20) of the order parameter 
equation with respect to Xh and taking into account Eqs. (3.17) and (3.29), one obtains 

A-KC) - ff§XAO + ^ 0. (4.48) 

where 


A XfiXmoi^Xji) T XfYiol^Xli) 1 2,XYno{xii') "h xfj 

= Xmo{Xh) - 2Xmo{Xh)X'f^{l/2\Xt,Xh). 


(4.49) 
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The derivation of the solution of the linear first-order ordinary differential equation 
(4.48) which meets the condition X^{f)\xt,Xh, Xmo) = 0 is straightforward and so we 
can express the local susceptibility dfx(C) through the order parameter Xm{C) fhe 
following explicit form 


.Y,(C) = A'i.(C) r XJfiffidw. (4.50) 

Jo [X'^{,w)] 

Starting from Eq. (4.45) one can also determine the function Xrnoixh)- Obviously, Eq. 
(4.45) can be written in the form 


^ O— X^o') (T/j, vY^^q) (4.51) 

OXh oXmO 

and since X^ {f)\xtiXhi X^q) = 0 one derives 


j^mo{Xh) 


Xfi^ XxnO^ 
dXmO Xfi, X^ffj 


(4.52) 


Thus, all the terms in Eq. (4.50) are completely determined only by means of the scaling 
function of the order parameter prohle X^ and its derivatives. Furthermore, Eq. (4.51) 
delivers an alternative analytical expression for Xy. provided one knows X^- 

The local susceptibility functions corresponding to the pair of order parameter 
prohles depicted in hgure 6 that are coexisting at capillary condensation curve are 
obtained using Eqs. (4.51) and (4.52) and are presented in hgure 7. From this hgure we 




Figure 7. The local susceptibility functions corresponding to the left and right order 
parameter profiles (see figure 6) coexisting at capillary condensation for Xh = —250 
and Xt = —33.8105. 


see that on the capillary condensation curve the local susceptibility might have either 
one ot two local maxima. The inspection of hgures 6 and 7 let us to conclude that in 
the case when the local susceptibility is characterized by two symmetrical local maxima 
they are centered, approximately, around the two gas-liquid interfaces in the system. 

4-2. Exact results for the local susceptibility profiles in the case of zero field 

When h = 0 one can determine the scaling function Xy^{z\xt) = X^{z\xt,Xh = 0) (see 
Eq. (4.44)) of the local susceptibility in an explicit analytical form [46]. One has 

W(Cki) = lAiKki) + c2i/'2(Cki). 


(4.53) 
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V'i(Ckt) = 




1 - 2dn 


.X 


m,0 


C;fc' 


and 


with 




X 




fc'(l - 2/c'^)E ( am (V fc' ] - ik'^Xm,oC 


k' 


+k'cn(t^C,-k' 


C2{Xt) = — 


k 

Ak'^k'^K (k) 


.X 


fc'2 + (i_2fc'2) an(t^-^C,k']2 


k'^K (k) + (F - A;'2) E (k) ■ 

Here E{k) is the complete elliptic integral of the second kind. 
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(4.54) 


(4.55) 




(4.56) 



Figure 8. Color online. Plot of the function Xj^(^|a;t) normalized per value of this 
function in the middle of the system for T = Tc- 

The behavior of the scaling fnnction for three different temperatnres - well 
below, at and well above is shown on figure 8. 


f.3. Exact results for the total susceptibility in the case of nonzero field 

According to Eq. (1.4), for the scaling function of the total susceptibility X {xt, Xh) one 
has 

X{xt,Xh)= [ X^{C\xt,Xh)dC, (4.57) 


wherefrom one obtains that 


X {xt,Xh) = 2 


[X^(C)-X^o][X^(C)-A] 

) [XUO? 


dC 


(4.58) 


The temperature behavior of the susceptibility for several fixed values of the field 
scaling variable are given in figures 9 and 10. Figure 9 illustrates the dependence of 
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X, 



-35 -30 -25 -20 


X, 


Figure 9. The scaling functions of the total susceptibility —236.0350005), i.e., 
near Tc^l (left), and X{xt,—237.24), i.e., near Tcap (right), normalized with the total 
susceptibilityX(0,0) at the bulk critical point. X {xt,Xh) shows a clear singularity at 
Tc^l- The small finite jump of X (xt,Xh) near Tcap happens when xt passes through 
the pre-capillary condensation line. 


the total susceptibility in the vicinity of the points and Tcap. The left sub-hgure 
there shows the variation of the total susceptibility in the vicinity of the critical point 
= (-25.6983391,-236.0350005) of the hnite system, while the right 
one shows its behavior around the capillary condensation point Tcap = = 

(—26.4025, —237.2395). Here is understood as the highest temperature at which 
the entire capillary hlls with liquid. Figure 10 represents the variation of the total 




X, 


X, 


Figure 10. The typical temperature behavior of the scaling function of the total 
susceptibility X{xt, —250) (left) and X{xt, —50) (right) well below and well above the 
critical value of the magnetic scaling variable x^f \ correspondingly, normalized with 
the total susceptibilityX(0,0) at the bulk critical point. 

susceptibility with the temperature for helds far away from The left hgure 

corresponds to Xh = —250 < x^fi^ and the hnite jump of X {xt,Xh) indicates the passing 
of Xt through the capillary condensation line. The curve on the right hgure demonstrates 
that X {xt,Xh), as expected, is smooth when the held Xh = —50 is well above x^f^\ 
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f.f. Exact results for the total susceptibility in the case of zero field 

Having in mind Eqs. (4.53) - (4.56), one can also determine [46] the scaling function 
X{xt) = X(a;t,0) of the total susceptibility x(xt), where 

X{xt) = Vl^X{xt), (4.59) 

with 1 / = 1/2 and 7 = 1 for the model considered. The corresponding result for X{xt) 

is 

, c.,(x,)IK(k) + K(k)-2E(k) 

=-iifvh-■ 

Here C 2 is given in Eq. (4.56) and Xt is to be determined from Eqs. (3.22) and (3.24). 



Figure 11. The scaling function of the total susceptibility X{xt) normalized per its 
value at the bulk critical point. 

The behavior of X{xt) is illustrated in hgure 11. One observes that X possesses a 
maximum above the critical temperature of the bulk system at Xt = 19.9678 and the 
value of the maximum is 1.09925 times higher than the value of the total susceptibility 
at the critical point. 

5. Discussion and concluding remarks 

In the current article we present exact analytical results for the temperature-held 
behavior of the order parameter prohle and for the behavior of the main response 
functions - local and total susceptibilities, for one of the basic and most studied models 
in the statistical mechanics - the mean-held Ginsburg-Landau model. We have 
studied the properties of this model under the so-called (-1-, -b) boundary conditions 
for a system with a him geometry in the case when both bounding the system surfaces 
strongly prefer the liquid phase of the conhned huid system. We studied both the critical 
regime tL^/'" = 0(1), hL^/'-' = 0(1), where t = {T — Tfi/Tc, v = 1/2, A/z/ = 3, and 
the capillary condensation regime T < Tc, h < 0. The basic new exact result is the one 
derived for the behavior of the scaling function of the order parameter prohle presented 
in Eq. (3.30), wherefrom one derives expressions (4.49) ~ (4.52) for the scaling function 
of the local susceptibility, and Eq. (4.58) for the total susceptibility. Analytical results 





Exact results for the temperature-field behavior of a mean-field model 


19 


for the behavior of these quantities were known before only for the h = 0 case. The 
behavior of the order parameter prohle for different values of T and h is visualized in 
hgures 1, 3, 5 and 6, of the local susceptibility - in hgures 7 and 8, and that one of the 
total susceptibility - in figures 9, 10 and 11. 

Based on the derived exact analytical expressions we obtained the phase diagram 
(see hgure 4, left) and the coordinates of the critical point Tc^l = = 

(—25.6983391, —236.0350005) which are in excellent agreement with those determined 
in [44], see Fig 13 therein [55]. We observed that along the coexistence line between Tc^l 
and the the capillary condensation point T^ap = = (—26.4025, —237.2395) 

the system jumps from a less dense gas to a more dense gas (see hgure 5) before switching 
on continuously into the usual jump from gas to liquid state in the middle of the system 
in the capillary condensation regime - see hgure 6. This is an unexpected theoretical 
result that calls for an experimental check-up — one shall see if it is a theoretical artifact 
of the considered model, or corresponds to experimentally observable phenomena. 

Closing this discussion, let us also mention that the mean-held solutions are exact 
for the critical behavior of systems with dimensionality d > 4 (apart from some 
logarithmic corrections for d = 4). Next, these solutions serve as a starting point 
for more sophisticated analytical techniques like the renormalization group calculations 
utilizing the e-expansion [30, 43, 50]. Thus, our results shall be also helpful for such 
future theoretical considerations. 
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Appendix. Some details on the numerical evaluations 

In the current appendix we present a numerical verihcation of Eq. (3.43) and provide 
some technical details related to the determination of the numerical results presented 
in some of the plots in the current article. 

As already stated in the main text, using Eq. (3.43) the difference AT between 
the energies of two states is determined by the difference between the two respective 
truncated energies up to terms of order 0(d), i.e. 

AT = ATi,(d)+0(d), (A.l) 

where 

AT S{xti Xhi Xmofi) T(xt, x/i, A^qi), (A.2 ) 

^h) Xm02i Tjj-(x^, X/i, , d). (A.3) 

The relation (A.l) means that i) d determines the precision with which we determine the 
energy difference AT between any two solutions of the oder parameter problem and ii) 
that hm^_,.o ATir(d) = AT. A numerical verihcation of Eq. (3.43) is presented in hgure 
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Al. This figure comprises two typical examples of the evolution of AStr as a function of 6 



logio(<5) logioW 

Figure Al. The evolution of the difference AEtr between the truncated energies of 
the two competing states for 10“^^ < S < 10“® at Xt = —80,a;?i = —331.89 (left) and 
Xt = —2Q,Xh = —236.55 (right). 

for Xt = —80, Xh = —331.89 (left) and Xt = —26, Xh = —236.55 (right). The choice of 6 is 
indeed important because choosing S too large can lead to a wrong conclusion about the 
energy difference of the competing states. As figure Al shows, the energy difference is 
negative for S = 10“®, but becomes positive at 5 = 10“^ and stabilizes below S = 10“®. 
Thus, one should carefully choose the truncation S in equations (3.42) and (3.43) in 
order to ensure that the difference of the truncated energies properly indicates which of 
the two competing states for a given {xt,Xh) combination is of less energy. To find such 
6 one, e.g., chooses 6 = 10“™ and computes the difference of the truncated energies for 
increasing integers m till this energy no longer changes sign. Figure Al depicts how the 
value of AStr approaches its limiting value AS with the decrease of S, in accordance 
with Eqs. (3.43) and (A.l). The insets therein show the variation of the difference 
AStr = logio — AStr(d)) with log]^Q((5). As we see, this dependence is linear 

which is exactly what it shall be expected on the basis of Eq. (3.43). 



-338 -332 -326 -236.61 -236.55 -236.5 


Xh 

■ gas-like; • coexistence; ♦ liquid-like 


Xh 

■ gas-like; • coexistence; ♦ denser gas-like 


Figure A2. The evolution of the difference AEtr between the energies of the two 
competing states for xt = —80 (left; the usual jump from gas-like state to liquid-like 
state), and Xt = —26 (right; a jump from a gas-like state to a denser gas-like state). 
The circles indicate the coexistence points for these temperatures. 
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At the end, let us present some clarification remarks about the procedure followed 
in the determination of the phase diagram given in figure 4. If, at any fixed {xt,Xh), the 
system possesses two competing states we determine the one which is stable by choosing 
that one with less energy. An illustration of the evolution of AStr with the change of the 
thermodynamic parameters governing the behavior of the system is shown in figure A2. 
The point {xt,Xh) for which AStr vanishes within the chosen precision 6 , is the point 
that belongs to the phase separation line of the phase diagram. The last implies that 
the phase diagram is, of course, also determined within that precision. In the current 
article, as stated above, we have worked with S = 10“^^. 
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